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Restricted equilibrium ensembles: Exact equation of state of a 
model glass 
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PACS 05 . 70 . Ce - Thermodynamic functions and equations of state 
PACS 64 . 60 . My - Metastable phases 
PACS 02 . 50 . Ga - Markov processes 

Abstract. - We investigate the thermodynamic properties of a toy model of glasses: a hard-core 
lattice gas with nearest neighbor interaction in one dimension. The time-evolution is Markovian, 
with nearest-neighbor and next-nearest neighbor hoppings, and the transition rates are assumed 
to satisfy detailed balance condition, but the system is non-ergodic below a glass temperature. 
Below this temperature, the system is in restricted thermal equilibrium, where both the number 
of sectors, and the number of accessible states within a sector grow exponentially with the size of 
the system. Using partition functions that sum only over dynamically accessible states within a 
sector, and then taking a quenched average over the sectors, we determine the exact equation of 
state of this system. 
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The glassy state is usually considered a non-equilibrium 
state of matter in the sense that the conventional 
Boltzmann-Gibbs treatment in terms of partition func- 
tions would not work, and yields very different proper- 
ties of the corresponding "equilibrium state" (e.g. of crys- 
talline quartz, and not of window glass). In this paper, our 
aim is to develop an equilibrium statistical mechanical de- 
scription of the glassy state using restricted phase-space 
ensembles. We consider a system constrained by dynamics 
to a restricted region of phase space. Within this region, 
it acts like an equilibrium system. This may be considered 
as an idealized description for metastable states such as 
supercooled liquids, or a glass. Focussing on the latter, we 
note that materials like window glass, for a given history 
of preparation, and at temperatures sufficiently below the 
glass temperature, have a well-defined macroscopic den- 
sity, velocity of sound, and specific heat. Then, over a 
time scale of microseconds to years, these materials are in 
some effective restricted thermodynamic equilibrium. 

In the restricted equilibrium ensemble corresponding to 
a glassy state, the partition function sum only extends over 



a restricted region of the phase space. As an illustrative 
example of these ensembles, we discuss a simple model. 
In particular, we determine the exact equation of state (a 
material -dependent relation between the density, pressure 
and temperature). To the best of our knowledge, this is 
the first time the exact equation of state for a model with 
short -range interactions and showing a glassy phase, is 
obtained. 

The model we study is a lattice gas in one dimen- 
sion, with a pair potential. The time evolution is as- 
sumed Markovian, with detailed balance, but in the glassy 
regime, the dynamics is non-ergodic, and phase space 
breaks up into a large number of disconnected sectors. 
We calculate averages of physical quantities in terms of 
partition functions within a sector, and then average this 
over sectors. 

The notion of restricted ensembles is not new. Many 
authors have discussed metastable states in terms of con- 
straints on the regions of phase space available to the sys- 
tem, i.e. phase space is broken into disjoint components, 
with crgodicity within components [TJ. The idea of com- 
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ponents, in the context of supercooled liquids and glasses 
was made more specific as inherent structures by Still- 
inger and Weber [2]. This gives the energy or free-energy 
landscape picture, which has been very useful in provid- 
ing a conceptual framework for the study of glasses. But 
the actual calculation of the partition function within a 
component, or the number of components is very difficult, 
and one usually just postulates the form of the distribu- 
tion function of these [3]. The study of specific models 
with Markovian evolution with detailed balance, show- 
ing break-down of phase space into disjoint sectors, was 
pioneered by Fredrickson and Andersen [3], and such ki- 
netically constrained models have been studied a lot re- 
cently [5]. However, usually the cases studied have a triv- 
ial Hamiltonian, and the main focus is on the relaxational 
dynamics. Our focus here is on the thermodynamics in 
the glassy phase. 

We should also mention earlier work on metastable 
states for systems with long-range couplings, as in the 
Sherrington-Kirkpatrick model of spin-glasses [6], or su- 
percooled liquids 0[5]. Theories like the mode-coupling 
theory of glasses mainly discuss the onset of glassiness 
from the liquid side [9]. An overview of recent work on 
supercooled liquids and structural glasses may be found 
in PHUT]. 

Our approach uses pico-canonical ensembles, discussed 
earlier in [12]. The glass transition is built into the sys- 
tem by assuming temperature-dependent rates of local dy- 
namic processes, some of which are set to be exactly zero in 
the glassy phase. We do not try to describe the approach 
to the glass transition, but focus on the near-equilibrium 
behavior away from the glass transition point, in the glassy 
phase. 

Definition of the model. — The system we study 
is a hard-core lattice gas on a semi-infinite line of sites, 
labelled by positive integers i, with i varying from 1 to oo. 
At each site i, there is an occupation number variable rii 
that takes values or 1, depending on whether the site is 
occupied or not. 

There are a total of N particles on the line. On the left, 
there is an immovable wall at i = 0. To the right of the 
rightmost ( the iV-th) particle, is a movable piston, whose 
position will be denoted by L, and which exerts a constant 
pressure p on the system. There is an attractive interac- 
tion between nearest neighbor occupied sites of strength 
J. For simplicity, we assume that the interaction between 
piston and particles of the gas, and between the left wall 
and the particles is the same as that between two particles. 
Then, with the convention no = Ul = 1, the Hamiltonian 
of the system is given by 

L-l 

H = - J njUj+i +pL (1) 

i=0 

We assume that the system evolves by Markovian dynam- 
ics, with the following rules: 
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Fig. 1: A schematic representation of the glassy and non-glassy 
regions of the control parameter space (p,T). The non-ergodic 
glassy regions are shown as shaded regions G. A denotes a 
path in this space, corresponding to particular preparation of 
the system. 



i) An occupied site with an empty neighbor can exchange 
position with the neighbor with a rate T\ exp(— /3AE/2). 
Here AE is the difference of energies between the final and 
initial configurations, and j3 is the inverse temperature. 
This is represented by the 'chemical' equation 01 ^ 10. 

ii) An occupied site with two empty neighbors on its 
left ( or right), can jump two spaces to the left ( right), 
with a rate r2 exp(— (3AE/2). This is represented by the 
equation 100 ^ 001. 

iii) The piston can exchange position with a neighboring 
empty site. This process is represented by the equation 
OP f± P0. The corresponding rate is T^exp{— (3AE/2). 

In our model, the externally controlled variables are the 
pressure, and temperature. The dynamics conserves the 
number of particles. Within the model, we are free to pre- 
scribe any functional dependence to the rates T\, T% and 
T3, and the long-time steady state of the system, and the 
thermodynamic properties do not depend on the precise 
values of these parameters. We postulate that F2 and T3 
are non-zero for all temperatures, and may be assumed to 
be independent of temperature and pressure, without any 
loss of generality. For T\(p, T), we assume that in two- 
dimensional control space (p, T) , there are regions where 
it is zero, and in the rest of space, it is not [fig. 1]. The 
regions where Fi(p, T) is zero will be called glassy regions. 

The dynamics of this model is somewhat similar to the 
model of diffusing reconstituting dimers (DRD model) dis- 
cussed earlier in [12] . In the DRD model, the only allowed 
transitions in the low-temperature phase are 110 f± 011, 
while in the high temperature phase, 01 ±=? 10 is also al- 
lowed. These moves are the same as the moves in the 
present case if the 0's and l's are interchanged. 

The DRD model was defined in the constant L ensem- 
ble, and the right wall could not move. It is not possible to 
define a satisfactory dynamics of a piston consistent with 
the dynamics of the DRD model. For example, consider 
the configuration 10001 1001 1000100P. Then, it is easily 
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seen that in the DRD dynamics, the piston cannot move to 
the left of the present position of the rightmost 1. Hence 
the DRD lattice gas is incompressible in the non-ergodic 
phase. 

Again, if we move the piston in the DRD model to the 
right by some macroscopic distance, only dimers can move 
into the new region. Thus, the long-time steady state is 
inhomogeneous, and shows a phase separation, with the 
left region at a higher density than the right, and pressure 
equilibration between different parts is not possible under 
the DRD dynamics. This makes the DRD model unsuit- 
able for studying the equation of state. The present model 
does not suffer from these problems. 

Sector decomposition and the restricted ensem- 
ble. — We note that each of the Markov rates here satis- 
fies the detailed balance condition, and the state with the 
probability of a configuration C is given by the standard 
Boltzmann -Gibbs weight 

Prob(C) = exp[-/3if (C)]/Z(p, T) (2) 

is a steady state. Here Z(j>, T) is the partition function of 
the system, given by 

Z(p,T) = J2eM-PH(C)} (3) 

c 

If T) is non-zero, any configuration of N particles 
can be reached from any other. Then, the measure given 
by Eq. ^ is the unique long-time steady state measure. 

If however, Ti(p, T) — 0, then the phase space breaks 
up into a large number of disconnected parts, called com- 
ponents, or sectors. There are 2 N disjoint sectors. This is 
easily seen: Define the position of the «-th particle., 
and Ai — Xi—Xi-i. Then the different can only change 
by ±2. We define m = A l {mod2). Then {r] t },i = 1 to N 
are conserved. Note that f?jv+i is not conserved. Clearly, 
the number of sectors is 2 N . We will denote a sector by S, 
and label these sectors by specifying the N binary integers 
{*}. 

The conditional probability of a configuration C in the 
long-time steady state, given that the system is in the 
sector S is given by 

Prob(C|S) = cxp[-/3i?(C)]/Z s (p, T). (4) 

The normalization factor Z§(p,T), called the pico- 
canonical partition function for the sector S, is defined 
as the sum of the standard Boltzmann weights, but only 
over configurations C that are in the sector §, 

Zs(p,T) = J2eM-PH(C)} (5) 

ces 

We consider the control parameters (p, T) being varied 
slowly along some curve A in the control parameter space 
(Fig. 1). The curve A may enter and leave the glassy 
regions several times. When the system is in a non-glassy 



region, there is a unique equilibrium state of the system, 
that is independent of the system history, and only de- 
pends on the current value of the control parameters. In 
the glassy regions, the restricted equilibrium state depends 
on (p, T), and also on the sector. If we vary (p, T), staying 
within the glassy region, the sector does not change. The 
relative probabilities of different sectors depend only on 
the point from which one entered the glassy region. 

A precise specification of the sector § requires N bits of 
information, and if N is large, this is neither possible, nor 
useful. In an experimental set up, S is specified in some 
general way in terms of how the system is prepared. We 
assume that in the beginning, the system is prepared in 
the ergodic region of the control-parameter space (p,T), 
and then brought to the desired state by moving along a 
specified path A adequately slowly. We will specify the 
macroscopic state of system in a glassy state by the last 
values of (p, T) in its history, when it was ergodic with 
Ti / 0, to be denoted by (p*,T*). Then, the probability 
that the system falls in the sector S is given by 

Prob(S) = Z s (p*,T*)/Z{p*,T*) (6) 

For any observable 0, the long-time average value 
within the sector § will be denoted by (O)g. This is given 
by 

(O} s = ^0(C)Prob(C|§) (7) 

ces 

Often, these can be expressed as appropriate derivatives 
of the free-energy, as in the standard equilibrium statisti- 
cal mechanics. For example the average length (£)§ of the 
system in the sector § is given by 

{L) s = -k B T^-\ogZ^{p,T). (8) 

The expected value of an observable (O) in an experiment, 
for a given history of preparation of the sample, is obtained 
by further averaging (0)§ over different sectors § 

(0)=^Prob(§)(0) s , (9) 

s 

where Prob(§) is the probability that the system will freeze 
into the sector §, for the given history of preparation. 

One can similarly calculate variance of the observable, 
or the sector-to-sector variation of the sector-mean (0}§. 
In the simple case we are considering, macroscopic quan- 
tities like the mean density can be shown to be self- 
averaging, and the variance of the sector-means are of 
0(N). The fluctuations, compared to the mean, are 
smaller by a factor TV -1 / 2 , for large N. If the relative 
fluctuations are not small for some observable, predicting 
its value in a particular experimental set-up, without any 
additional information about the experiment, is clearly 
not possible. In such cases, one can only determine the 
probability distribution of such an observable. 

In particular, the mean free energy in the glass phase is 
defined as the free energy, averaged over sectors. We are 
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working in the constant (p, T)-ensemble, and the appro- 
priate free energy is the Gibbs free energy $(p, T). The 
average value of <&{p, T) in the glass phase, averaged over 
sectors, will be denoted by ($g(.P; T)). This is defined as 

(*G(p,T)> - -fc B T^Prob(§)logZ s (p,T) (10) 



Note that ($ G (p, T)) implicitly depends also on (p*,T*). 

Equation of state. — We now determine the equation 
of state for this system. Consider first the ergodic case 
when Ti ^ 0. In this case, any allowed configuration can 
be reached from any other under the Markovian dynamics 
of the system. 

We define xn+i to be the position of the piston L. In 
terms of the variables {A^}, the Hamiltonian can be writ- 
ten as 



N+l 

H = V[-J<5 A ,,i +pA ?; 



z=l 



(11) 



Clearly, in this case, the different Aj are independent ran- 
dom variables, we get the partition function of the system 
Z(/3,p) in the constant temperature and pressure ensem- 
ble as [13] 



Z(p,T) = w(u,x) 



N+l 



(12) 



where we have used the notation u = e^ J , and x — e 
and 

„2 



w(u, x) 



\ u 5&1 x A = UX + 
^ 1 - X 

A=l 



J- 



(13) 



The mean spacing between the molecules (A), which is 
the inverse of the density of the lattice gas, is given by 



(A) = x— log w(u, x) 



w(u, x) 



2x-x 2 



forT^O, (14) 



and the mean energy per particle (E(u,x)} is given by 



(E(u, x)) = —Ju 



1-x 



, for Ti ^ 0. (15) 



We can similarly calculate Zg(p, T), for any given sector 
8. Define 

u 5a - 1 x a (16) 



Wodd 



A odd 



and a similar equation for w even , where the sum is only 
over even values of A. Then we have 



Wodd = ux + 



1-x 2 



and 



Weven = X 2 / (1 - X 2 ) 



(17) 
(18) 



We denote the number of zero rji's by r(S). It is easily 
seen that 

„ r (S) AT-r(S)/ . s n s 

Z S = W e y J n W odd '(Wodd + Weven) (19) 

Under quenching, rji takes the value with probabilities 
w even /w = a* independent of i, and the value 1 with 
probability 1 — a* . 



(20) 



X* + x* 2 + u*(l — x* 2 ) 



Then, the probability of a sector S with exactly j odd rji's, 
is given by 

Prob(S) = a*\l -a*) N ~ j (21) 
The variable r(S) is distributed as a binomial distribution 

Prob(r(S) = j) = N Cj a* 3 {I - a*) N ' j (22) 

Thus the distribution of r(S) is sharply peaked, with max- 
imum at a*N, with a width of order 0{N 1 / 2 ). 

A straight forward calculation gives the mean spacing 
between the particles in the low-temperature phase as 

d 

(A) = x— [a* loguw„(w, x) + (1 - a* ) log w od d{u,x)]. 

(23) 

Using the expressions for w dd and w eV en (Eqs. 12 and 
13), we finally get 



(A) = (1 - a*) 



-, 2x 2 (i-x 2 y 



u+x^l-x 2 )- 1 



for Ti = 0. (24) 



This equation, along with Eq. (1141) gives the mean vol- 
ume per particle as a function of u and x, and hence is 
the equation of state of the material. It depends on the 
history of the system through the parameter a*, which 
depends on onset of the glassy state. From Eg. (|2"2l . the 
sector-to-sector variation of r(E>)/N is of C(iV -1 ' 2 ), and 
hence fluctuations in (A)§ are also small for large systems. 

Similarly, we can calculate the mean energy per particle 
(E(u, x)) by taking derivatives with respect to u, and we 
get 



(E(u,x)) 



~ J "L a *L ^ r ^°- ( 25 ) 



u + x 2 j (1 — x 2 ) ' 



As an example, consider the simple case J = 0. In this 
case, there is no energy, and the density is a function of 
only one variable Pp. In this case, it is easily seen that the 
equilibrium density is given by 



Peq = 1 - e 



(26) 



We assume that the high-density phase is not ergodic, 
say whenever p > 1 — x* . Then, it is easy to see from 
Eq. ([2~4"|) . that in the non-ergodic phase, the dependence is 



Pglas 



(l-e- 2 ^) [2a* + (l-a*)(l + e- 2ft5 )] 1 (27) 
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Fig. 2: Plot of density p as a function of the scaled pressure 
Pp for J — 0, when the regime /3p* > log(2) is assumed to be 
non-ergodic. The point of onset of non-erdodicity is denoted 
by G in the figure. 

where a* — x* /(l + x*). 

We have plotted p eq and p g i aS s in Fig. [5] for x* = 1/2. 
We note that the density is continuous at p — p* , and the 
glass is less compressible than the corresponding equilib- 
rium state. 

Entropy. — The statistical mechanical entropy of a 
macroscopic system in equilibrium, i.e. one described by 
the Gibbs ensemble, when the state state is discrete, and 
the Gibbs ensemble assigns probability pi to the ith state, 
is given by 

" (28) 



S = Pi log ft 



This entropy coincides, to leading order in the size of the 
system, with the thermodynamic entropy of Clausius, and 
is the same for different ensembles, for ordinary systems 
( we are ignoring gravitational interactions). The set of 
microstates used is always a restricted one, in practice. 
For example, when considering argon gas below 10 3 °K, we 
do not consider the ionized states, or the excited nuclear 
states. There is some controversy in the glass literature 
[T^] about the appropriate choice of the set of microstates 
over which one should sum in Eq. (|28"|) . In the following, 
we shall assume that the sum extends over all states which 
could be reached, given the thermal history of the system. 

There are L C N for each L, and L = N + 1, N + 2, 

This entropy can be expressed as a sum of two terms. 
Firstly, we do not know which sector the system has fallen 
into. The corresponding entropy, usually called the frozen 
entropy of the system, and is given by 



S 



frozen 



= - 5^ Prob(§) log Prob(S) 

s 



(29) 



This is easily computed in our case using Eq. (|21j) . and we 
get frozen entropy per particle Sf rozen , in the thermody- 
namic limit of large N given by 



For x* = 1/2, we have a* — 1/3, and the frozen entropy 
per particle is log 3 — flog 2 « 0.636, which is a large 
fraction of the total entropy per particle at the glass point 
s* = log(4) 1.386. 

Note that s frozen only depends on (u*,x*), and not on 
u or x explicitly. Hence in taking derivatives with respect 
to u or x, it docs not contribute, and the mean energy 
or pressure are same as would be computed from pico- 
canonical partition function. 

The second contribution to entropy comes from the 
many possible microstates within one sector. This de- 
pends on the sector, and its calculation involves the pico- 
canonical ensemble. This will be denoted by S P i CO , and is 
given by 

S P ico(S) = -J2 Prob(C|§) logProb(C|§) (31) 

This can be expressed in terms of the picocanonical par- 
tition function, for large N, to leading order in N, as 



Jpico \ 



log Z S -/?|^ log Z S (32) 



Let fj is the conditional probability that a particular 
A takes the value j, given whether j is even or odd. 
Clearly, f x = ux/w odd , f 2 j+i = x 2j+1 /w odd for j > 1, 
and f 2j = x 2 i- 2 (l - x 2 ). Using the fact that Prob(C|§) is 
a product of probabilities of different Aj's, it is straight- 
forward to calculate the per site entropy, s P i COl averaged 
over the sectors S, and we get 



j even 



(l -a*) 



j odd 



(33) 

The controversy in literature about the entropy of 
glasses, in our model, relates to the question whether to in- 
clude Sf roze n in the thermodynamic entropy. This seems 
to us a matter of taste. One can argue that if we start 
from the high temperature phase, and cool down along 
some annealing path A (Fig. 1), then the microstates 
in different sectors should contribute to entropy. On the 
other hand, if we prepare the system in some sector S, 
then other states not consistent with this initial prepara- 
tion would not contribute. It is clear is that Sf rozen will 
not change if we change control parameters staying within 
the glassy phase. This could be a change of temperature 
in the glassy phase, and whether we add Sf rozen or not, 
we would get the same specific heat in the glassy phase. 

We note finally that the 1-d lattice-gas model with near- 
est neighbor couplings can be seen as a renewal process, 
i.e. the probability that separation between the ith par- 
ticle and (i + l)th particle being j is independent of i. 
Denote this by Cj. A straightforward calculation gives 



Sfrozen = ~G* loga* - (1 - a*) log(l 



(30) 



fjCt*, for j even; 
fj(l — a*), for j odd. 



(34) 
(35) 
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With this distribution, the total entropy per site is given 
by 

s = Sf rozen + Spi CO = ^ ^ Cj logCj. (36) 

3 

Discussion. — One important drawback of this par- 
ticular model deserve mention. It does not have a crys- 
talline phase at low temperatures, and consequently there 
is also no supercooled liquid phase. In the model, there 
is a direct transition from the ergodic liquid to the non- 
ergodic glass phase. 

Of course, Ti becoming exactly zero in the glassy 
regime is an idealization, used to make the problem more 
tractable. It is not qualitatively different from other ide- 
alizations, such as infinite heat reservoirs with weak cou- 
plings etc. that are routinely used in equilibrium statisti- 
cal mechanics. 

One can ask how does the relaxation rate between dif- 
ferent sectors depend on T\, when Ti is very small. In 
this case, if one looks at the process as a change in the 
sector- labels {rji}, with time, the process is like a sym- 
metric exclusion process in this space, where 01 f± 10 with 
rate Ti. From general arguments, the relaxation time of 
this process should grow as L 2 /T\. 

In general, for a non-equilibrium system, with a Hamil- 
tonian -ff(A), in contact with a heat bath at a fixed tem- 
perature T, if a control parameter in the Hamiltonian is 
changed from A, by a small amount to A + AA, one can 
write AE, the change in the mean energy E of the system 
as 

AE = AJ2 Pi^i = [PiAEi + EiApi] . (37) 

i i 

Here the first term may be identified as the work done on 
the system AW, and the second as the amount of heat 
added AQ [15]. In our model, if we go across the glass 
transition, neither pi, nor Ei undergo any change. Hence, 
there is no work done on the system, or heat exchanged 
with reservoir, as we cross the ergodic-glass boundary. 

It is quite straightforward to generalize this model. For 
example, instead of step-size 2, we can consider step size 
3, or 4 etc.. We can also consider longer-ranged interac- 
tions. This does not change the sector decomposition, but 
calculation of partition function within a specified sector 
becomes more complicated, but still reducible to a finite 
matrix diagonalization, if the range of interaction is finite. 

We can also construct models in higher dimensions. 
Consider, for example, a lattice gas on a square lattice, 
with a pair-wise additive interaction potential of finite 
range. At high temperatures, we assume that each particle 
can diffuse to its nearest neighbor at a rate satisfying the 
detailed balance condition. Now, we break up the lattice 
into cells of size k x k, ( k is any integer) and at low tem- 
peratures, a particle is only allowed to hop within its own 
cell. Then, in the low temperature phase, the number of 
particles in each cell gets frozen in time, and the number 
of such sectors grows exponentially with the volume of the 



system. However, calculating the pico-canonical partition 
functions in such models has not been possible so far. 

If there is a sequence of ergodicity breaking transitions, 
(say corresponding to a smaller region G\ inside G, in 
Fig. 1), then the state will depend on more parameters 
a*, a|, . . ., and would show a much more complicated de- 
pendence on the history of the system. 
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